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Abstract 

Previous work has shown that static multi-black hole solutions of 
higher-dimensional Einstein-Maxwell theory do not possess smooth hori- 
zons. We show that the lack of smoothness is worse than previously 
demonstrated. We consider solutions describing multiple black holes on 
a common axis. In five dimensions, the metric is generically twice, but 
not three times, continuously differentiable at the horizon. The Maxwell 
field is generically continuous, but not differentiable, at the horizon. In 
more than five dimensions, the metric is once, but not twice, continuously 
differentiable, and there is a parallely-propagated curvature singularity 
at the horizon. The Maxwell field strength is again continuous, but not 
differentiable, at the horizon. 

1 Introduction 

It has long been known that four-dimensional Einstein-Maxwell theory admits 
static multi-black hole solutions. The existence of these solutions can be at- 
tributed to a balance of electromagnetic and gravitational forces. The (degen- 
erate) horizon of each black hole in these solutions is known to be analytic 

ID- 

String theory motivates the study of multi-black hole solutions of higher- 
dimensional Einstein-Maxwell theory. Perhaps surprisingly, it turns out the 
horizons of such solutions are generically not smooth. This was first suggested 
by Gibbons et al [5] and later proved by Welch [3] for the case in which the 
black holes lie on a common axis. Welch considered a timelike geodesic along 
the axis as it approaches the horizon of one of the black holes. He calculated 
the components of the Riemann tensor in a parallcly propagated frame. For 
a two-black hole solution, he showed that second derivatives (with respect to 
proper time) of some of these components diverge at the horizon. This proves 
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that the Riemann tensor is not C 2 (twice continuously differentiable) at the 
horizon and hence the metric is not C 4 at the horizon. Below we shall present 
a much shorter proof of this result. 

Stating things more precisely, the work of Welch proves that there does not 
exist a C 4 extension of the metric through the horizon. However, the question 
still remains: how smooth is the solution? To answer this question, it is nec- 
essary to construct a coordinate chart that covers the horizon. One such chart 
was constructed in 0]. The resulting metric is C° but not C 1 . Hence previous 
work has demonstrated that there exists a C° extension, but not a C 4 extension. 
The aim of the present paper is to sharpen these results. 

To our surprise, it turns out that higher-dimensional multi-black hole so- 
lutions are even less smooth than proved by Welch. By considering non-axial 
geodesies, we shall show that the metric at the horizon of a two-black hole 
solution is not C 3 in five dimensions. We shall then introduce Gaussian null 
coordinates to construct a coordinate chart in which the metric is C 2 in a neigh- 
bourhood of the horizon. Hence the metric is C 2 , but not C 3 . The Maxwell 
field strength is C°, but not C 1 , at the horizon. 

In higher dimensions, things are even worse: the metric is C 1 , but not C 2 , 
in more than five dimensions. In fact, we shall show that there is a parallely 
propagated curvature singularity at the horizon. This sounds problematic: one 
usually demands that the metric should be C 2 in order that the Einstein equa- 
tion makes sense. However, the Maxwell tensor is still C° and hence the RHS 
of the Einstein equation is C°, so the Einstein tensor must be C° even though 
there is a curvature singularity at the horizon. Therefore the Einstein equation 
is still satisfied at the horizon. 

What lies behind the horizon? When dealing with a non-smooth horizon, 
this question need not admit a unique answer. We shall show that there are in- 
finitely many different ways of matching the solution across the horizon onto an 
interior solution, whilst preserving the degree of differentiability just discussed. 
This phenomenon also occurs for four-dimensional cosmological multi-black hole 
solutions [5]. 

As in the cosmological case [5] , the degree of differentiability at the horizon 
of a given black hole can be increased by introducing extra black holes in an 
appropriate (finely-tuned) manner [3] . One can even make the horizon of one 
black hole analytic. However, the horizons of the other black holes will still be 
non-smooth unless one considers an infinite array of black holes. Surprisingly, 
this "smoothness enhancement" is only possible in d = 5: there are no d > 5 
solutions with a metric smoother than C 1 at the horizon. 

This paper is organized as follows. In section [2 we describe the general 
multi-black hole solution with which we will be working and present a brief 
argument that demonstrates that the horizon of a two-black hole solution is not 
C 4 . In section [3] we analyze d = 5 multi-black hole solutions in detail. The 
d > 5 solutions are discussed in section 2) Some technical details are given in 
the Appendices. 
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2 Multi-centre co-axial solutions 



The solutions we are considering in this paper are those of ci-dimensional Einstein- 
Maxwell theory, with action: 



/ 



d^xV=g[R-[J r -^r)F IH> F^) (I) 



The static multi-centre solutions to this theory have the following metric: 

ds 2 = -H- 2 dt 2 + H 2 '^ (dr 2 + r 2 d6 2 + r 2 sin 2 8dn 2 _ 3 ) (2) 

where dQ 2 ,_ 3 denotes the metric on the unit (d — 3)-sphere. The Maxwell field 
strength is F = dA, where 

A = -H- X dt. (3) 

The function H is harmonic on R d ~ 1 , with poles corresponding to the locations 
of the event horizons of the black holes. We have chosen our coordinate system 
in the metric above such that the harmonic function for the multi-centre co-axial 
solution depends only on r and 8. Thus the multi-centre harmonic function for 
a line of N black holes in our case is: 

N 

TT _ 1 , Mi , \p M» 

r d ~ 3 ^ (r 2 +af- 2a,r cos 9W~ 3 )/ 2 { ' 

where the a» denotes the position along this axis of the i-th black hole, and /U, is 
a mass parameter. We shall investigate smoothness of the horizon at r = 0. To 
this end, it is convenient to expand H in (axisymmetric) spherical harmonics: 



H= + £ h n r n Y n (cos 6), (5) 



Mi 

7 "t 

where the harmonics ^(cos^) are given by Gegenbauer polynomials: 



r„(cos(9) = Cn 2 \cos6) (6) 



with the first few being 

d-3 

C 2 =1 

(7^ = (d-3) cos (7) 
c ^ = (^3) ((rf _ 1)cos 2^_ 1) 

The coefficients h n are given by 

^ = ^o + EuTC- (8) 



i=2 



For generality, we shall take h n to be independent parameters until we need to 
specify to particular examples. 
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2.1 Short argument for non-smoothness 

The argument of [3] involved calculating derivatives of components of the Rie- 
mann tensor in a parallely propagated frame. Here we shall present a much 
shorter argument that allows one to reach the same conclusion, namely that the 
metric is not C 4 at the horizon of a two-black hole solution. 

We start from the assumption that the metric admits a C k extension through 
the horizon of one of the holes, with k > 2, and that this extension admits a C 1 
Killing vector field V which coincides with d/dt outside the horizonQ Note that 
C 1 is the minimal assumption for Killing's equation to make sense. However, 
we can now show@ that V must actually be C k . As V is a Killing vector, we 
have the following relation: 

v Q .v b y c = -Ri a V d (9) 

where R^ ca is the Riemann tensor. As the metric is C k , the Riemann tensor 
must be C k ~ 2 , so at least C° as k > 2. Hence the RHS of © is C° and therefore 
so must be the LHS, which implies that V must be C 2 . Plugging this back into 
©, the right hand side must now be C min ^ 2 ' k ~ 2 \ so by the previous argument 
the Killing vector must be (7 mm ( 2 > fc - 2 )+2 = £*nun(4,k)_ Repeating this argument 
we will eventually learn that V is C mm (™> fc ) where n > k, so by induction we 
arrive at the conclusion that V must be C k . 
Now consider the norm of V: 

V^V = -H- 2 (10) 

Since both the metric and V are C k , it follows that H~ 2 must be C k at the 
horizon. We can use this to determine an upper bound on k by considering the 
differentiability of H~ 2 along axial null geodesies. 

Instead of the spherical polar coordinates of ©, we can use cylindrical polar 
coordinates (p, z, fids) so that the metric has the form 

ds 2 = ~H- 2 dt 2 + H 2 ^ d -V [dp 2 + dz 2 + p 2 dn 2 _ 3 ) . (11) 

For a two-centre solution we have 

TT _ -I J Ml | ^2 (v2 s 

(p 2 + z 2 )( d - 3 )/ 2 [p 2 + {z- a) 2 )( d " 3 )/ 2 ■ V ; 

Consider a future-directed null geodesic approaching the origin along the posi- 
tive z-axis. Such a geodesic has p = 0. From energy conservation we have (for 
appropriately normalized affine parameter A) 

l = " 3 < 13 > 



An alternative set of assumptions for the d = 5 case is motivated by supersymmetry. 
The solution i ]2|.">[ i admits a super-covariantly constant spinor e, from which one can construct 
a Killing vector field V a ~ £7a£ [5] which coincides with d/dt. We could assume that the 
extension through the horizon preserves supersymmetry, i.e., e extends through the horizon. 
If the Maxwell field is C° (and the metric assumed C k , k > 2 as in the main text) then e will 
be C 1 and hence so will V . 

2 We thank P. Chrusciel for this argument. 
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and from the null condition we have 



dz d-t |=| / z d - 3 z d - 3 rn \ d - 3 

dX \ /ii //i(a - z) d A J 

Solving this equation for small z (with z — > 0+ as A — * 0—) gives 

z{\) = (rf-3) 1 ^- 3 VP^(-A) 1 /( d - 3 ) 

+ (id - 2)(d - 3)^ M 7 1/(d - 3)2 (l + a 3 -V 2 )(-A)^ 

+ ^-^(d - i)^^^a 2 - d {-\)^ + 0((-X)^) (15) 
2d — 5 

This implies that, along this geodesic we have 

H- 2 = (d- 3)V7 2/(d - 3) A 2 + (d - 3) 2 (d - 2) Ai - 2/(d - 3) (l + a 3 -V 2 )A 3 

- 2(d " 1 ^_" 5 3) ^^(-A)^ + 0((-A)^F) (16) 

This is an expansion in powers of A 1 /^ -3 ^. For d = 4, such an expansion is 
analytic in A, as expected from pQ. However, if d > 4 then the third term is not 
C 4 at A = 0. Hence H~ 2 is not C 4 at the horizon so we must have k < 4, i.e., 
the metric is not C 4 , which was the conclusion of Welch 0. In later sections, 
we shall consider the behaviour of invariants along non- axial geodesies and show 
that this leads to the conclusion that the metric is even less smooth. 

2.2 Gaussian null coordinates 

In later sections we will employ Gaussian null coordinates to extend the metric 
through the horizon. In this section we give a brief description of how such 
coordinates are constructed. More details can be found in [7]- We shall restrict 
attention to extensions that preserve the R x SO(d— 3) symmetry of the metric 
1(2} . An argument of the previous section shows that the associated Killing fields 
will have the same degree of smoothness as the metric. 

Consider a single component of the event horizon TLq . The intersection Ho 
of this with a spatial hypersurface is topologically a sphere S d ~ 2 . Introduce 
coordinates x % on S d ~ 2 . Let V denote the generator of time translations (V — 
d/dt in the coordinate of ([2])). V is tangent to the null geodesic generators 
of Hq (this follows from the fact that V generates a symmetry and is null on 
the horizon). Define a coordinate v on Ti^ to be the parameter-distance from 
Ho along integral curves of V. This defines a coordinate chart (v,x l ) on a 
neighbourhood of Ho in ■ 

Now let U be the unique (past-directed) null vector field satisfying U -V = 1 
and U ■ d/dx l = on TCq. Let j(v,x l ) be the null geodesic that starts at 
the point with coordinates (v, x l ) in Hq and has tangent U there. Gaussian 
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null coordinates are defined in a neighbourhood of Hq by ascribing coordinates 
(v, X,x l ) to the point affine parameter distance A from Hq along 7(1;, x l ). The 
metric in these coordinates takes the form 

ds 2 = -H(X,x)~ 2 dv 2 + 2dvdX + 2Xh i (X,x)dvdx 1 + h ij (X,x)dx i dx j . (17) 

Note that V — d/dv in these coordinates, and hence g vv is determined by 

y2 = _ H -2 

We can be more explicit about the coordinates x % on S d ~ 2 . It is convenient to 
take x % = (9, Qd-3), where (6, Od-3) are the limiting values of (9, Qd-3) along 7 
as 7 approaches Tig". The SO(d—2) symmetry implies that angular momentum 
is conserved along 7. Hence, since we have initial condition dfld-3/dX = at 
A = 0, we must have dQd-3/dX = 0, i.e., Qd-3 is constant along 7. Hence 
fld-3 = However, 9 and 9 will not agree. In these coordinates, we will 

have 

ds 2 = -H-' 2 dv 2 +2dvdX+Xf 1 (X, 0)dvdO+f 2 (X, 9)d9 2 + J ff 2/(d_3 V 2 sin 2 9dfl d - 3 

(18) 

where now r = r(A, 9) and 9 = 9(X, 9) are determined once we know the 
geodesic 7. 

The problem of bringing the metric to the Gaussian null form (|18p amounts 
to solving the geodesic equation to determine ^(v, 9, ^-3) using the original 
coordinate system of ([2]). As we have already seen, Qd-3 is constant along 7. 
Energy conservation, together with U ■ V — 1 gives 

§—*>,«). (19) 

The null condition then reduces to 

= -H 2 + H 2 '^ (f 2 + r 2 9 2 ) , (20) 

where a dot denotes a derivative with respect to A. The final equation we need 
is the geodesic equation for r: 

r-H^d r H+-^—H- 1 r 2 d r H — -^—H- 1 r 2 9 2 d r H~r9 2 + -^—H- 1 r9d e H = 0. 
a— 3 a— 3 a— 3 

(21) 

The problem reduces to solving (f2T)]) and (f2"Tj) simultaneously to determine 
r(A, 9) and 9(X, 9), where the initial conditions are 

r(o,e) = o, fl(o,e) = e, ^) =o. (22) 

Gnce r and are known as functions of A and 9, equation (|19[) can be integrated: 
t = v-T{X,Q), T(A,9) = /" H(X,Q) 2 dX, (23) 
with u arising as the constant of integration. 
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2.3 Single black hole 

To illustrate all of this, let us first examine the single centre solution. The 
harmonic function in the metric of ^ is now simply 

ff W = l + ^3" (24) 

In this case, the solution is spherically symmetric. This additional symmetry 
implies that 9 = along 7. The constraint (|2"0)) reduces to 



which can be integrated to give 

d-3 



1 ■ (<i 3) 

( M 1/(d - 3) +A) -„ 



(26) 



Expanding out for small r we see that 

l/(d-3) 



Equation (j23j) is now 



r~ ((d-3) M (d - 4)/(d - 3) A) . (27) 

Writing iJ as a function of A gives 

# = — A ■ 28 

( M V(<i-3) + A) d - 3 - M 

i = v - J H(X) 2 dX (29) 
and hence 

dt = dv- H(X) 2 dX. (30) 
Finally we obtain the metric in the Gaussian null coordinate system: 

ds 2 = -H~ 2 dv 2 + 2dvdX + ( M V(d-3) + a) 2 (d9 2 + sin 2 6d0^_ 3 ) . (31) 

The metric components are analytic functions of A so this defines an analytic 
extension of the solution through the horizon at A = to negative values of 
A. In this case, spherical symmetry implies that g v Q = (and the Gaussian 
null coordinates coincide with Eddington-Finkclstcin coordinates). However, 
the multi-centre solution will take the more general form ([T? 
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2.4 Strategy for multi-black holes 



To find the functions r(A, 0) and 8(X, 0) for the multi-centre solution we sub- 
stitute the power series expansions 

r=J2 C m(®> h n)X n/{d - 3) 

m=1 (32) 

m—1 

into equations (f2"0)) and (|2"Tj) . The motivation for the choice of fractional powers 
in these expansions comes from the leading order behaviour of r as a function 
of A for a single black hole, given in equation (|2~T|) . The coefficients in these 
expansions are determined using computer algebra. 



3 Multi-centre solution in five dimensions 
3.1 Determining the geodesic 

We begin by considering the d = 5 solution, for which we can write the har- 
monics in the form 

sin ((n + 1)0) 

F n (cos6>) = — . (33) 

Upon solving the geodesic equation with the expansions (j3"2"|) as described in the 
previous section, and imposing the intial conditions (|22[) . we find the functions 
r(A, 0) and 8(X, 0). The first few terms are 

r = V2dW + — ^A 3 / 2 + 0(A 2 ), 

r ( 34 ) 

e = e _2_V2h^ne x3/2 + 0{x2y 

Mi 

More details, and further terms in the expansion, are given in Appendix A. The 
function T(A, 0) of equation (f2"3")) is also given there. 



3.2 The metric is not C 3 

To investigate the smoothness of the metric, we consider the area of the 2-sphcre 
orbits of the SO (3) symmetry in the geometry, which is a scalar invariant of the 
solution. The metric on these S 2 's is obtained by restricting the full metric to 
the space spanned by the SO (3) Killing fields. Since the Killing fields have the 
same differentiability as the full metric, it follows that the metric on these S 2 's 
must also have the same differentiability as the full metric. The area of S 2 is 

A 2 = Hr 2 sin 2 6. (35) 
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We can determine how A 2 varies along the null geodesic 7 as it approaches the 
horizon by substituting our expansions for r and 9 along the geodesic into this 
expression. This gives 



A 2 = jUi sin 2 8 + sin 2 9A + (hi - 4/j 2 /«i sin 2 9) sin 2 9A 2 



_ i^/^fWecosOA^ + C^A 3 ). (36) 



The presence of a term of 0(A 5 / 2 ) indicates that this quantity is not C 3 at 
A = 0, and we infer that the metric cannot be C 3 at the horizon unless /13 
vanishes. However, the explicit expression for h n JS} shows that /13 ^ for a 
2-centre solution. If there are more than 2 centres then /13 is still non- vanishing 
unless the parameters are finely tuned to make /13 vanish. We conclude that 
generic multi-centre solutions do not admit C 3 horizons. 

3.3 Transforming the metric 

We now want to construct a coordinate system in a neighbourhood of the horizon 
in which the metric is C 2 . The first step is to convert the metric §2$ outside 
the horizon to Gaussian null coordinates as explained above. The fact that the 
metric is not C 3 at the horizon will be reflected in the presence ofO(A 5 / 2 ) terms 
in the metric. 

In practice, r(A, 9) and 9(X, 9) are expressed as infinite expansions in A. In 
order to be as explicit as possible about the coordinate transformation, we shall 
truncate these expansions, keeping just as many terms as are required to make 
the metric coefficients C 2 functions of A. This amounts to keeping just the terms 
written out explicitly in Appendix A. Of course, by performing this truncation 
we are no longer dealing with exact Gaussian null coordinates as defined above, 
but with " nearly Gaussian null coordinates" , in which the metric differs from 
the Gaussian null form (fT5|) by terms of order A 5 / 2 . For example, there will be 
non- vanishing components g\\ and g\Q of order A 5 / 2 . 

The multi-centre d = 5 metric in these coordinates in a neighbourhood of 
the black hole horizon at r = is 



+ (fit sin 2 6 + 2VmT^o sin 2 9 A + (hi - 4/i 2 Mi sin 2 9) sin 2 OA 2 + C(A 5/2 )J dn 2 2 



The metric is not well-defined at 9 = 0, it. However, there is no indication in 
the above expansion that these are anything other than the usual coordinate 
singularities of spherical polar coordinates. 




+ 



+ 0(A 5/2 )dA 2 + 0(X 5/2 )d\dQ + 0{X 5 / 2 )dvd@ 
(txi + 2^h a X + (h 2 Q + 8h 2 m sin 2 9) A 2 + 0(A 5/2 )) d® 2 



(37) 



9 



3.4 Extending through the horizon 

We will extend the solution through the horizon as follows. We assume that the 
interior metric takes the form with the same coordinate ranges (e.g. r > 0) 
but with a different harmonic function H (r, 9): 

ds 2 = -H- 2 dt 2 + H (dr 2 + r 2 d6 2 + r 2 sin 2 Od^l), (38) 

where H is chosen to agree with H at leading order: 

oo 

H=^ + J2~ h nr n Y n (cos6). (39) 

n=0 

Recall that we are assuming that the Killing field V that generates time trans- 
lations can be extended through the horizon. We choose the time-orientation of 
the interior solution so that V = d/dt. 

We assume that the interior Maxwell field is as in ([3]) but with a sign change 
(which obviously still gives a solution): 

A = H-i-dt. (40) 

We will convert the interior metric to the "nearly Gaussian null" form, and then 
attempt to match it onto the nearly Gaussian null form of the exterior metric 
given above. This will impose some restrictions on the coefficients h n . 

In the exterior region, A > and A is the affine parameter along a past- 
directed geodesic 7. It is convenient to define a parameter A in the interior 
region to be the affine parameter along a geodesic 7 that obeys all the same 
conditions as 7 except that 7 is /irfure-directed (which amounts to replacing 
the condition U ■ V = 1 with U ■ V = — 1 above). We then construct Gaussian 
null coordinates as before. The only difference is that we will have g \ = — 1 
whereas g v \ = +1. To determine 7 we proceed as before. The equations are 
exactly the same except for a sign change in (fl~9"l) because the geodesic is now 
future directed. This does not affect the equations governing r and 6 hence 
the expansions for r and 9 are exactly the same as before but with A — + A and 
h n — * h n . Hence we can define nearly Gaussian null coordinates as before, with 
the coordinate v defined by 

t = v+f(\,e), r(A,e) = J H(X,e) 2 dX. (41) 

We can determine T(A, 0) by replacing A — > A and h n — > h n in T(A, 0). The 
metric in these coordinates is obtained from the exterior metric (|3T[) by the 
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same replacements and the sign change g v \ = 1 — > g vX = — 1: 

ds 2 = ( -—A 2 + l^A 3 + 0(\ 7/2 ) | dv 2 + (-2 + 0(X 9/2 ))dXdv 
V ^ Mi / 

+ 0(A 5/2 )dA 2 + C(A 5/2 )dAd9 + C(A 5/2 )dwde 

/tii + 2V7IT^oA + (/i 2 , + 8/i2Mi sin 2 g) A 2 + 0(A 5/2 )) d9 2 

(fix sin 2 9 + sin 2 9 A + (M - 4h 2 fii sin 2 0\ sin 2 9A 2 + 0(A 5/2 )) dflj 

(42) 



Finally, we define the coordinate A for the interior region by A = —A. It is then 
clear that the interior metric (|4"2")) for A < matches onto to the exterior metric 
(|37[) for A > in a C 2 manner (i.e. up to C(|A| 5/2 )) provided that 

ho = —ho h 2 = h 2 - (43) 

The other multipole moments h n are unconstrained, hence infinitely many inte- 
rior solutions can be matched onto a given exterior solution so that the metric 
is C 2 at the horizon. Somewhat surprisingly, h\ is unconstrained. 

3.5 The Maxwell field 

We may now check the degree of differentiability of the Maxwell field. In the 
exterior region, the potential in the nearly Gaussian-null coordinates is 

A = -H- l dt = -H- l {dv - d x Td\ - d e TdG). (44) 

Expanded out for small A, this is 

a = (--La + ^a 2 + 32 ^c° s e A5/2 + 3 \ dv 

\ Mi S^ 7 J 

+ (V2 x - 1 + 3ho + h^^ hl cos QA i/2 + 0{x) ) dX 
\ Mi 4 5 J 

_ ^2^1/4^ s . n QA 3/ 2 + 0(A aj j de (45) 

The singular C(A _1 ) term and the constant term in A\ are obviously pure 
gauge. In the interior region we have 

A = H- 1 dt = H- 1 (dv + d k fd\ + d e fdQ) (46) 
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which gives 



A = 




+ 




A 5 / 2 + C(A 3 )^ dv 



Upon defining A = — A for the interior region it is clear that the Maxwell field 
strength is C° in these coordinates. However, it is not C 1 because F\@ — 
C( | A | ) unless the solution is finely tuned so that h\ — and we choose 
h\ = (and ho = —h as in (|43)) ). in which case the Maxwell field strength 
becomes C 1 . 

3.6 Smoothness enhancement 

We have seen that the metric of a d = 5 multi-black hole solution is generically 
C 2 but not C 3 at the horizon, and the Maxwell field strength is generically C° 
but not C 1 . However, we have commented that the differentiability of the metric 
can be increased if hs = and that of the Maxwell field increases if hi = 0. 
This increase in differentiability for certain configurations was also discussed 
in [3] for the particular case of a three-black hole solution, where the masses 
and locations of the black holes can be chosen in such a way that the central 
black hole horizon may be made more differentiable. A similar effect had been 
observed previously for cosmological multi-black hole solutions in d = 4 [5] • 

In general, the lack of smoothness arises from the fact that r ~ A 1 ' 2 , so odd 
powers of r are not smooth functions of A. Hence the degree of differentiability 
increases if we arrange for the lowest odd multipole moments in the expansion of 
H to vanish. If we we arrange the black holes so that all odd multipole moments 
vanish then the horizon (of the black hole at r — 0) becomes smooth, in fact 
analytic. We shall now discuss this case in more detail. 

Assume that the sources are arranged in such a way that there is a reflection 
symmetry in the plane 8 = ir/2. In other words we have N = 2M + 1 black 
holes with parameters (/j,i, 0), (/Zj, ±aj), i — 2 . . . M + 1. Then ft-2n+i = for all 
n. We write H has 



H = H + H, H = ^ + h a , ff = ^/i 2n r 2n y 2n (cos0). (48) 



We can now construct an analytic extension through r = using essentially the 
same method as in d = 4 [T]. Define A > by (compare (|25p) 



oo 



71=1 








1/2 



(49) 
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which gives 




1/2 



(50) 



Define v by (compare (|30|) ) 



dt = dv- H (XfdX. 



(51) 



The solution in coordinates (v, A, 9, f^) is then analytic in A at A = and can 
therefore be extended through the horizon to negative values of A. In the interior 
of the black hole we can define a new coordinate r as 



and a new coordinate t by (|5Tj) to bring the interior solution to the form (|38I40[) 
with a new harmonic function 



Note that this obviously obeys the conditions f4"3")) required for the extension 
to be C 2 , but now we see that analyticity uniquely determines all of the higher 



This solution is analytic at the horizon of the black hole at r = 0. However, 
it will not be smooth at the horizons of the other black holes. It is interesting 
to ask whether one can construct a multi-black hole solution for which all black 
holes have analytic horizons. The above analysis, reveals that this will be the 
case if the sources arc reflection-symmetric about every black hole on the axis. 
It is easy to see that this implies that there must be infinitely many black 
holes present, equally spaced on the axis of symmetry, with masses alternating 
between two values \x\ and fi2- Such a solution with /ii = /12 was studied in 
[4], where a periodic identification was imposed in order to obtain a solution 
describing a single black hole localized on a Kaluza-Klein circle. If \i\ 7^ ^2 then 
one can periodically identify in order to obtain a solution describing two black 
holes localized at antipodal points of a Kaluza-Klein circle. (Obviously one can 
also construct solutions for which the two black holes are localized anywhere 
on the circle but the horizons will be analytic only when the black holes are at 
antipodal points.) 

4 Multi-centre solutions in d > 5 

Having warmed up with the case d = 5, we now proceed to look at the smooth- 
ness for d > 5. Using the power series ansatz given in section [2 we can solve the 
geodesic equation order-by-order as before, with d allowed to take any integer 
value greater than 5. Upon substitution into the geodesic equation (|21[) , with 
the null constraint (|2H)l and the initial conditions (f2"2"|) imposed as before we find 




(52) 



H(r,0) 



H(ir, 9). 



(53) 



coefficients h 2n too. 



c 2 = c 3 = ... = c d _ 3 = 
bi = b 2 = ■ ■ ■ = b d - 3 = 



(54) 



13 



and the first few non-zero terms in the expansion are 



(d-4), 



c d -2 = ^-^(d - 3)^Mi ( " 3) ho 
c ^-i = — £ — — W ^1 cos6 I 55 ) 

2d — 5 



d-2 



b d -2 = ~{d- 3)5-5 Mi C5=^ ftl sin0 

= ^(d-3)^//p^/i 2 sin(20). 

With these we can check the smoothness of the area of the S d ~ 3 along the 
geodesic. The area is 

A d _ 3 = H^r 2 sin 2 6 (56) 
which expands out for small A to give 

2 1 3d— 11 

A d -3 = Mr 3 sin 2 9 + 2/z<^/i sin 2 6 A + /Ltp^ (d - 3)^U 2 sin 4 eA( 1+ ^) 

+ 0(A( 1+ ^)), (57) 

2d-6 

where we have neglected terms of order 0(X d - 3 ) as these are only of the same 
order as A( 1+2 /( d ~ 3 » for d = 5. 

For d = 5 the term of order A 1+2 /( d ~ 3 ) is smooth so the first term containing 
a non-integer power of A is of order As, which implies that the metric is not 
C 3 as we saw above. However, if d > 5 then the term of order A 1+2 /( d ~ 3 ) does 
spoil smoothness, and implies that the metric of a d > 5 multi-centre solution 
is not C 2 at the horizon. Note that it is not possible to arrange the black holes 
so that the coefficient hi vanishes (equation {8J) hence there is no analogue of 
the "smoothness enhancement" that is possible in d = 5. 

Using the same procedure as that used for d — 5, we can find the expansions 
(|32|) for d > 5. With these it is a straightforward task to introduce Gaussian 
null coordinates as described above. The solution in these coordinates is given 
in Appendix B. The extension through the horizon is performed in exactly the 
same manner as for d = 5 (section 13.41) . Again we assume that the interior 
solution is described by a new harmonic function H (r, 9) and Maxwell potential 
given by 



A = H-t-dt 



n=o (58) 



and then we determine the restrictions on h by matching this interior solution 
(expressed in Gaussian null coordinates) with the exterior one given in Appendix 
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B. The analysis goes precisely as described in section [5^1 We have already seen 
that a C 2 extension of the metric does not exist, so the best we can obtain is a 
C 1 extension, which imposes the condition 

h = -h Q . (59) 

The higher order coefficients hi and /12 appear in the metric above order A and 
so are not fixed by the requirement of matching the solutions in a C 1 manner. 
The Maxwell field strength in these coordinates is C° but not C 1 . 

The fact that the metric is not C 2 suggests that there may be a curvature 
singularity at the horizon. To see that this is indeed the case, note that some 
components of the Riemann tensor in the Gaussian-null coordinates diverge at 
the horizon. For example. 

R xee = 2 ^Agee)gee - (coffee) 2 ^ 
4gee 

which expands out as 



R xeev = (d - l)(d - 3) 2 /( d - 3 VP^/i 2 sin 2 eA^ + 0(a£S). (61) 

This diverges on the horizon for d > 5. By considering the Riemann tensor in 
a parallely propagated frame, we can show that this is a genuine curvature sin- 
gularity. We construct frame vectors parallely propagated along a null geodesic 
with tangent vector U = d/d\ as follows. Set e° = U. Then, on the horizon, 
define e 1 = V and e 2 = g^^d/dQ so that e° • e 1 = 1, e° • e 2 = e 1 • e 2 = 0, 
(e ) 2 = (e 1 ) 2 = 0, (e 2 ) 2 = 1. Extend e 1,2 off the horizon by demanding that 
they are parallely propagated along the geodesic: U ■ We 1,2 — 0. This preserves 
the orthogonality relations. Then, in this basis, we have 

R 0221 = R abcd e 0a e 2b e 2c e ld = g^Rxeev + ■■■ (62) 

where the ellipsis denotes subleading termsH From this we see that R 0221 di- 
verges at the horizon so there is a parallely-propagated curvature singularity 
there. 



5 Discussion 

We should mention two potential loopholes in our work. First, we have con- 
sidered multi-black hole solutions for which the black holes lie on a common 
axis. The exterior solution has R x SO(d — 3) isometry group, and we have 
considered only extensions that preserve this symmetry. It is conceivable that 

3 Here we have used the fact that the leading non-smooth terms in the metric arc 
C)(A( d - 1 )/( d - 3 )) and hence the leading non-smooth terms in the Riemann tensor in Gaus- 
sian null coordinates are 0{\^~ d ^^ d ~^ ). Therefore no component of the Riemann tensor is 
more divergent than R\qq v . 
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smoother extensions exist, but break some of the symmetry. However, for d > 5, 
this cannot be the case because we have shown that the solution has a curva- 
ture singularity at the horizon, and hence the metric cannot be C 2 there for 
any extension. For d — 5, we expect that a similar argument could be used 
to prove that some components of the curvature tensor in a parallely propa- 
gated frame must have divergent derivatives at the horizon, proving that no C 
extension of the metric exists. This would require constructing the parallely 
propagated frame beyond leading order. However, a simpler argument shows 
that the Maxwell tensor cannot be extended in a C 1 manner: the basis compo- 
nent F 02 = F Xa e 2a = Fxee 26 = 0{X 1/2 ) is not C 1 along the geodesic (here we 
used e 2v = from U ■ e 2 = 0). 

Second, our extension through a horizon makes use of spherical polar co- 
ordinates, and therefore exhibits singularities at the poles, where the horizon 
intersects the axis of symmetry. There is no indication that these are anything 
other than coordinate singularities (the axis itself is regular between the black 
holes) but a more rigorous treatment would require investigation of this issue. 
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A Expansions for r and 9 for d = 5 

The full expansions for r(A, 0) and 8(X, 0) are now given, up to the order 
required to eliminate from the metric all divergent terms and fractional powers 
up to A 5 / 2 . The coefficients in these expansions were obtained by substituting 
the ansatz ([521) into equations Q2Up . ([21]) and solving order by order. Note that 
the initial condition dd/dX = implies that 62 — 0. It is convenient to relax 
this condition, and take the limit 62 — > at the end of the calculation. This 
is because the equations exhibit degeneracy when 62 = 0. For example, 63 is 
determined by plugging the expansions into (|20j) and examining the coefficient 
of A 1 / 2 , which gives an equation of the form 62(^3 — . . .) = 0. Hence if 62 7^ 
then this determines 63. However, if 62 = then one has to go to higher order 
to determine 63 (with the same result). 
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r = V/ 2 + ka ■ n \*/*+l hl cos 9A 2 + 32 ^ ° SC 6 ^ = 3fe " A 5//2 

5 48V2^ /4 



4(20/i 3 /ii cos(29) - /i /ii) cos 8 3 
+ 35Vm 



25ftg + 64^i((86cos(26) - 89) hi + 25(1 + 2cos(29) + 2cos(40))/i 4 Mi) A 7/ 2 

1600V2/Lii /4 

(3/ig/ii + 24h h 3 fj, 1 cos(20) + 28 / u i ((4 + 8 cos(40))h 5 /zi - 50/n/i 2 sin 2 9)) cos 9 ^ 4 

63/ii 

+ ^ (33600(1 + 2cos(29))/i^2Mi - 7875/1 4 , 

1612800V2M1 7 

+ 2304/i Mi((6683 - 6610 cos(29))/i 2 + 245(1 + 2 cos(29) + 2 cos(49))h 4 A*i) 
+ 1024 A1 2 ((140cos(29)-7455+7630cos(49))^+72((194cos(49)-116cos(29)-58)/ii/i 3 
+35(1 + 2 cos(29) + 2 cos(49) + 2 cos(60))/i 6 /xi))) A 9/2 (63) 

2\/2/ii sin 9 3/2 . , 0/ ~ na2 

= 9 A' 372 - 3h 2 sm(29)A 2 

Mi 

+ V ?., (17ftofti sin 9 - 8h 3 ^i sin 9 - 24h 3 Hi sin(39)) A 5/2 

2 

(5/i 4 Mi sin(29) + 10/i 4 Mi sin(49) - 9hj sin(29) - 5/i ^ 2 sin(29)) A 3 



5VMi 

(64) 

T = -^A- 1 + - W/HlnA + fc + ^h^ 4 cosQX 1 ^ 2 
4 4 v 5 

+ C HM + -h 2 m csc9sin(39)^ A+— ^ /4 (19h /ii+40/i 3 Mi cos(29)) cos9A 3/2 
\ 16 3 / 35 

+ 777^-7= (25/ig + 75O(l + 2cos(20))/i o /i2A«i +4/zi((469- 106 cos(29))h 2 
400»//Zi 

+175(1 + 2 cos(29) + 2 cos(49))/i 4 Mi)) A 2 + tjt ((345ftgfti 

315VW 

+4272/i /i 3 Mi cos(29) + 224//i((40 - 13 cos(20))/nft 2 + 8(1 + 2 cos(49))h 5 Mi)) cos 0) A 5/2 
' (71400(1 + 2cos(20))/i^2Mi - 1575/iq + 144/i Mi((943 + 866 cos(20))/i 2 



100800^1 

+2205(1 + 2 cos(20) + 2 cos(49))/i 4y ui) + 256/i 2 (70(51 + 32 cos(20) - 11 cos(49))/i 2 

+9((416 + 832 cos(20) - 228 cos(49))^i^ 3 + 105(1 + 2 cos(29) + 2 cos(49) + 2 cos(60))/i 6 A«i))) A 3 

(65) 

with k an arbitrary constant. 
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B Gaussian null coordinates for d > 5 

The metric for d > 5 in the Gaussian null coordinate system is given by 

j2 / (d-3) 2 (d-2){d-3) 2 h ,o 2(d-l)(d-3) 1 ^ i , . , ™-s „-,,., 8--t A , s 

V (2d-5)/4 /(d - 3)2 V 



+ Uf" 3 +2 A t 1 d - 3 /i A- (d-3) — [i[ d - 3) h 2 sin 2 e\— +0{\—i) ) de 

2 1 3d- 11 

Hf* sin 2 6 + h sin 2 OA - /j< d - 3) " (d -S)$=*h 2 sin 4 6A^ +C(A^) ) df7 d _ 3 . 



(66) 



For simplicity we have assumed that the expansions for r and have been taken 
to infinity, allowing complete elimination of the g\\ and g\Q components. The 
Maxwell potential is 

(d - 1/(d - 3) A + {2 - d f- 3) ^ d -Vho\ 2 
l-d 



2d -5 



(d-3) — fi{ d - 3) ^cosGX— +0{\— ) ) dv 

d-2 



+ ((d-3)-Vl /(£i " 3) A- 1 + 



2(d-3) 



ho 



d-1 
2d -5 



(d - 3)^=3 Mf -33 /iicosOA— +C(A— ) dA 



d-1 
2d - 5 



(d-3)^>^- 3) "/iisin9A^ + 0(A^) ) dQ. (67) 



As in d = 5, the singular and constant terms in A\ may be removed with a 
gauge transformation. 
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